Method, instrument and system for performing verification measurements and corrections of magnetic fields in magnets for nuclear magnetic resonance imaging machines

ABSTRACT

The invention is directed to a method, an instrument and a system for performing verification measurements and corrections of magnetic fields in magnets for Nuclear Magnetic Resonance imaging machines. According to the invention, test elements are used whereof Nuclear Magnetic Resonance images are detected. Then, these actually detected images are compared with the theoretical images which would be obtained in ideal conditions. The number, magnitude and position of correcting magnetic charges are determined from the differences between the actual image and the theoretical image. The invention is also directed to a specific test element and to a Nuclear Magnetic Resonance imaging machine.

[0001] The invention is directed to a method for performing verification measurements and corrections of magnetic fields in magnets for Nuclear Magnetic Resonance imaging machines, comprising the following steps:

[0002] making a test member having such characteristics as to generate predetermined and known Nuclear Magnetic Resonance images;

[0003] obtaining a Magnetic Resonance image from said test member;

[0004] comparing the known theoretical image with the detected image and/or the received signals wherefrom the detected image is reconstructed with the corresponding set of signals, related to the known theoretical image in any of their processing or imaging steps;

[0005] determining the differences between the detected image and the known image or between the set of detected signals and the set of signals corresponding to the known image;

[0006] determining correction parameters, i.e. the correcting magnetic charge and/or the number of correcting magnetic charges having a predetermined value, or the volumes of ferromagnetic material and their position on the magnetic structure.

[0007] Nuclear Magnetic Resonance imaging machines use electromagnetic echoes of previously excited nuclei to retrieve therefrom information for imaging. In order to obtain echo signals allowing to reconstruct images sufficiently corresponding to reality, at least a considerable number of the nuclei of the material being examined must be oriented with substantially parallel spins. To this end, static magnetic fields are used which must be comparatively intense and constant within a predetermined volume corresponding to a body part to be examined.

[0008] At the same time, additional magnetic fields are applied, the so-called gradients, which are used to select certain sections of the body part to be examined and to create a parameter for identifying the signal received from individual portions of the section under examination so that signals can be ordered for two- or three-dimensional imaging.

[0009] The tolerances required for a good correspondence between reality and the reconstructed image are very small, being of the order of a few millionths deviations from nominal values. Moreover, since the magnetic fields in use are relatively intense, the constructions of magnets have a considerable size, whereby the physically huge structure hinders the required construction accuracy. This is particularly relevant for the so-called permanent magnets, in which the magnetic field is generated by permanently magnetized materials and not by induction generators.

[0010] After manufacture, magnets are submitted to fine calibration of the magnetic field, the so-called shimming, which is designed to correct construction inaccuracies of the individual parts of the magnet, and includes measuring the static field by a probe, detecting aberrations with respect to nominal values and disposing magnetized or ferromagnetic correcting elements in several appropriate areas of the magnet structure as determined from the differences between the nominal values and the actual values of the field within the volume designed to accommodate the body or the body part under examination. Shimming is performed during fabrication at the manufacturer's site. This technique is described in greater detail, for instance, in patent application SV98A000015 filed by the applicant hereof.

[0011] Nevertheless, frequently and unavoidably, upon installation at the customer's site, manufacture settings appear to be changed or perturbed. This is especially caused by stresses acting on the big metal masses which form the magnetic structure, especially as this structure has a modular construction, or made of several parts fastened together.

[0012] Further, the personnel using the machines is not sufficiently qualified to verify their operating accurateness during use, except in a few exceptional cases. Hence, little aberrations might remain unnoticed.

[0013] Due to this, a problem arises to make functional testing of machines easier, at least below a certain accuracy degree, in an easy and safe manner, even for personnel less qualified than the personnel charged of fabricating and testing the machines at the manufacturer's site.

[0014] Document EP-A-0230027 discloses a method and a device according the above described method and which attempts to solve the above disclosed problems. Also U.S. Pat No. 5,545,995 and U.S. Pat. No. 5,055,791, disclose method of the kind described above.

[0015] The known disclosed methods, suggest solutions to the problem of testing and correcting imaging aberrations and anomalies in which the image of the phantom is determined by a combination of several coefficients of the theoretical mathematical function which describes the field and which combination of coefficients includes coefficients of low and high order. Carrying out the correction of the filed under these circumstances requests very long calculations which are not necessary to achieve a practical level of precision of the field. Indeed in the disclosed solutions the structure of the phantom is not suited to the structure of the function describing the filed being normally a spherical polynomial expansion. The images of the phantom and the aberrations thereof from the theoretical image are not determined only by a small number of low order coefficients of the polynomial expansion. Thus as stated above the mathematical comparison of the real image of the phantom and of the theoretical image and the calculation of the correction to be applied to the filed are very complicated and time consuming.

[0016] The invention has the main object to allow fast functional testing of the magnetic structure of Nuclear Magnetic Resonance machines, particularly after installation, in a relatively fast and simple manner.

[0017] The invention has the further object to allow an at least partly automatic correction of aberrations and anomalies occurring after installation without requiring special laboratory equipment, which is complex and costly and may be only used by very highly knowledgeable and specialized personnel.

[0018] The invention also aims at setting the bases for implementation of several automation degrees of the correction to be performed, with the smallest number of operation steps which might be carried out by an average specialized technician, specifically trained in the maintenance of one or more specific machines.

[0019] The invention achieves the above purposes with a method for performing verification measurements and corrections of magnetic fields in magnets for Nuclear Magnetic Resonance imaging machines, in which

[0020] A mathematical theoretical function describing the field is choosen, particularly a polynomial expansion, preferably a spherical polynomial expansion,

[0021] A phantom is provided with a structure suited to the said function describing the filed and which produces one or more selected images, each of these images being correlated to one or to a limited number of selected low order coefficients of the function describing the filed, such that the differences between a theoretical image of the phantom and the real image of the phantom depends only to the one or to the limited number of selected coefficients of the function describing the magnetic field.

[0022] Thus by selecting a specific mathematical structure for describing the field, for instance and preferably but without limitation, a representation of the field by spherical harmonics, the above method may be refined so that several orders and degrees of harmonics, and hence of field coefficients may be examined individually.

[0023] In this case, the method provides that the test member has elements which do not emit Nuclear Magnetic Resonance signals and are defined transparent in the following text and claims , which are related to a definite harmonic or to precise coefficients of a certain order of a mathematical description of the field by a polynomial expansion, preferably by spherical harmonics, there being provided means for application of so-called reading gradients of the magnetic field, which only detect the echo signal along certain directions, said directions being selected in such a manner as to suppress the contributions from the magnetic field described by coefficients other than the ones being examined.

[0024] Preferably but without limitation, the describing polynomial expansion is selected in spherical harmonics. Nevertheless, other polynomial expansions may be more or less suitable for the purpose, also depending on the symmetries of the field and hence of the magnetic structure.

[0025] In accordance with a preferred embodiment of the method, the latter includes the following steps:

[0026] a) detecting the image of the predetermined test element;

[0027] b) symmetrizing deviations between the actually detected image and the nominal, ideal one, with respect to the center of the image and/or to the predetermined origin of the coordinate system for the mathematical description of the magnetic field;

[0028] c) defining a curve of symmetrized deviations having the relevant field coefficients as a variable;

[0029] d) determining a polynomial for approximating the curve of symmetrized deviations, deriving from the differences betweem the actual image and the ideal, theoretical image;

[0030] e) determining the coefficients based on the system of the two consistent equations of steps c) and d);

[0031] f) computing, from the mathematical description in spherical harmonics of the field, such number, magnitude and position on the magnetic structure of the correcting elements as magnetic charges or volumes of ferromagnetic material, as to bring the measured values of the field coefficients back within the nominal values corresponding to the homogeneity characteristics of the field required for detecting useful Nuclear Magnetic Resonance images;

[0032] g) manually positioning such charges on the magnetic structure.

[0033] Possibly, once correction has been performed, steps a) to f) may be repeated to verify the effectiveness of said correction and, when necessary, an additional correction step g) may be performed.

[0034] The elements for verifying the different coefficients may be progressively mounted on the test member or fixedly provided thereon.

[0035] Each of these elements corresponds to at least one harmonic or to at least one set of coefficients and are made consistently with the selected mathematical description of the field and/or eventually to the symmetries of the magnetic structure.

[0036] According to a first embodiment, the method provides a manual graphic comparison between the detected images and the theoretical images wherefrom the mathematical computation is performed to determine correction data.

[0037] In this case, the machine will have a theoretical comparison image, for instance stored therein.

[0038] Alternatively, the method provides that the machine determines quantitatively the differences between the ideal image and the detected image by computation of the distance in pixels of the point or portion or area of the detected image from the position of the corresponding point or area of the ideal image. In this case, the difference may be easily changed to an aberration or deviation quantitative value with respect to the theoretical image and directly used in the machine thanks to the processing hardware contained therein and to software for effecting comparisons, tests and error diagnostics. Hence, computation software, also stored in the processing hardware of the machine may compute the required corrections and indicate correction absolute data.

[0039] In these conditions, the operator would only have to physically position the correction charges, based on the magnitude suggested by the machine, in the position/s automatically determined by the machine.

[0040] In a specific embodiment, the method provides the use of a test element whose non-detectable portion has the form of a rib, stem or baffle whose longitudinal orientation is parallel to the static magnetic field, i.e. disposed in the axial section plane of the magnetic field when the latter is described with spherical harmonics. In a description by a polynomial expansion with spherical harmonics, this stem is related to (2 0) and (4 0) cosine coefficients.

[0041] From two crossed baffles which form an angle of 109.472°, symmetrically with respect to one of the two axes which describe the axial section plane of the magnetic field, the coefficients of the (2 2) cosine and of the (2 1) sine are obtained.

[0042] Similarly, from two crossed baffles which form an angle of 90°, symmetrically with respect to an axis which describes the coronal section plane of the static magnetic field the coefficients of the (2 2) sine are obtained, whereas from the same baffle configuration, with baffles disposed in the sagittal section plane of the static magnetic field, the coefficients of the (2 1) cosine are obtained.

[0043] When these coefficients are verified and, if needed, corrected, they are sufficient to ensure a good accuracy of the field to obtain the tolerances required for detecting useful Nuclear Magnetic Resonance images.

[0044] It has to be noted that (x y) typically defines the indexes of the polynomial expansion.

[0045] The invention is also directed to a test element including, separately or in combination, the above baffles.

[0046] In the preferred embodiment, such test element includes a central rectilinear baffle which is crossed, with reference to the plane parallel to a peripheral edge, by two crossed baffles, passing through the median area of the rectilinear baffle and forming, symmetrically with respect to the rectilinear baffle, an angle of 109.472°, whereas two additional crossed baffles are provided in the plane perpendicular to the rectilinear baffle and containing the other sides thereof, which form, symmetrically with respect to said rectilinear baffle, an angle of 90° and intersect the rectilinear baffle in the central area.

[0047] Advantageously, the baffles consist of walls and are made of plastic, particularly of the so-called Plexiglas.

[0048] In accordance with the illustrated preferred embodiment, the different baffles made of a material which emits no Nuclear Magnetic Resonance signal are mounted in a fixed manner inside a plastic container filled with a liquid which emits a Nuclear Magnetic Resonance signal, e.g. water or other liquids or solutions.

[0049] It has to be noted how the creation of contact surfaces between materials differing with respect to Nuclear Magnetic Resonance signals is relevant to generate image contrasting imaging areas, having predetermined shapes. Obviously, said preferred construction allows the method to be implemented even when the materials composing the baffles and the filling are inverted.

[0050] The box or container have a filler opening with sealing means.

[0051] The advantages of the present invention are apparent from the above disclosure. Thanks to the verification and correction method, the homogeneity accuracy of the magnetic field may be periodically controlled in an easy manner, manually and/or with different automation degrees. This feature is primarily important for obtaining Nuclear Magnetic Resonance images with a good correspondence with reality, hence it is the basis to obtain diagnostically usable images.

[0052] Verification and, when needed, correction activities, may be automated to such an extent that they may be performed by an average specialized operator, without requiring any intervention of highly specialized personnel, which involve high costs both for the manufacturer and for the customer.

[0053] Further, verification and correction times are drastically reduced, thereby limiting downtime, which causes discomfort for users.

[0054] The required process and devices are relatively inexpensive, mainly involving the addition of software procedures in existing machine elements, which only need to be enhanced, and the provision of a set of correction permanent magnets, as well as of a phantom. The latter is the only construction element which has to be specially fabricated for an easy application of the method.

[0055] However, the phantom itself is in no way necessary, the verification and correction procedures being performable, by following the above method steps, with hand-made elements or improvised members.

[0056] Further improvements of the invention will form the subject of the subclaims.

BRIEF DESCRIPTION OF THE DRAWINGS

[0057] The characteristics of the invention and the advantages derived therefrom will appear more clearly from the following description of a non limiting embodiment, illustrated in the annexed drawings, in which:

[0058]FIGS. 1A and 1B to 4A and 4B show schematically by a tridimensional view and a plan view on the specific extension plane, the test elements for verification of specific field coefficients according to a spherical polynomial expansion.

[0059]FIG. 5 shows a flow chart of the method according to the invention.

[0060]FIG. 6 shows the ratio of two coefficients as a function of the angle formed by the test baffle in the specific plane.

[0061]FIGS. 7 and 8 are two views along two axes, of the composite test element according to a preferred embodiment of the invention.

DETAIILED DESCRIPTION OF THE DRAWINGS

[0062] With reference to FIGS. 1 to 5, the method of the invention provides Nuclear Magnetic Resonance imaging of a test element and comparison of the deviations of said real detected image with a nominal, ideal image of the test element. In FIG. 5, the process is shown with reference to the test element example of FIGS. 1A and 1B.

[0063] In FIG. 5, I denotes the ideal image and I′ the actually detected image of the test element 1, according to FIGS. 1A and 1B. I+I′ denotes the comparison between the two images. This is obtained by superposing the two images I, I′. In fact, verification might as well be carried out visually, and aberrations be detected by manual measurement on the monitor. However, depending on the desired automation degree, such graphically represented verification may be obtained automatically by pure computation, whereas the visual image may only be a graphic representation of the information obtained from computations.

[0064] Then, aberration values, i.e. substantially the differences between points P1 and P1′ to Pn and Pn′ are used for computing the corrections to be applied to the magnetic structure 1.

[0065] Such corrections are provided in the form of magnetic elements M1, M2, Mn, having a unit value or being a fraction or a multiple of said unit value, or in the form of one or more elements with a predetermined volume of ferromagnetic material, which have to be positioned at predetermined locations of the magnetic structure to correct any aberration of the magnetic field, here of the static magnetic field denoted by Bz in the figures.

[0066] Preferably, the magnetic field is described with the help of a polynomial expansion in spherical harmonics which relate coefficients of different orders to combinations of sine and cosine functions. This description allows to produce test elements whose image contributions only relate to defined coefficients describing the magnetic field in spherical harmonics and have very simple geometries and constructions for this purpose.

[0067] Thus, for instance, the test element of FIGS. 1A and 1B allows separate detection of (2 0) cosine and of (4 0) cosine coefficients.

[0068] The following description will refer to the different geometries shown in FIGS. 1A and 1B to 4A and 4, while explaining the process in greater detail.

[0069] Referring to FIGS. 1A and 1B, the test element consists of a rectilinear baffle 4, which is positioned in the plane parallel to the magnetic field Bz, obviously between the two poles 2, 3 for generating the magnetic field.

[0070] The rectilinear baffle is oriented along the Z axis. In this condition, the tesseral coefficients with a non-zero degree give a null contribution on the axes.

[0071] The ideal image which should be obtained is denoted by I in FIG. 1. If the field is inhomogeneous, the actual image might be the one denoted by I′ in FIG. 5.

[0072] In this case, by adding the distortion of each point along an axis perpendicular to the field and selected in accordance with the direction of the reading gradient, for instance the x axis, to the point symmetric thereto with respect to the origin 0, i.e. the central point of the baffle, the contributions from the other symmetries are suppressed (see FIG. 5 I+I′, IS). By assigning this sum divided by two to each point and to its symmetric counterpart (P, P′, P1″, P′″, P2, P2′, P2″, P2′″), a curve (resulting from symmetrized construction) is obtained in a chart being a function of the distance from the center, which describes the distortion related to the relevant symmetry. This curve consists of powers of r, where r is the distance from the center of the image, i.e. from the origin, in coordinates normalized to the radius of a sphere having a predetermined diameter, e.g. 80 mm.

[0073] The inhomogeneities of the magnetic field Bz in the points symmetrized on the y axis are described by the following function: δ = a₂₀^(cos)r² + a₄₀^(cos)r⁴ + O(r⁶)

[0074] The distortion of the image is proportional to inhomogeneity according to the following formula, indicated as (1): $\begin{matrix} {{\Delta \quad {X\lbrack l\rbrack}} = {\frac{\Delta \quad {B\lbrack{MHz}\rbrack}}{f_{samp}\lbrack{MHz}\rbrack}{{FOV}\lbrack l\rbrack}}} & (1) \end{matrix}$

[0075] Where l corresponds to the unit of image distortion normalized to the radius of the above defined sphere.

[0076] Considering a number of 256 pixels, the above formula gives: ${\Delta \quad {X\lbrack{pixels}\rbrack}} \cong {{2.56 \cdot 10^{- 4}}\frac{f_{magnet}\lbrack{MHz}\rbrack}{f_{sampling}\lbrack{MHz}\rbrack}\left( {{a_{20}^{\cos}r^{2}} + {a_{40}^{\cos}r^{4}}} \right)}$

[0077] By approximating the distortion of symmetrized points as a function of their distance from the center of the magnet by a polynomial function like the following:

y=β ₂₀ _(^(r)) ²+β₄₀ _(^(r)) ⁴

[0078] The desired coefficients are obtained, i.e. (2 0) cosine and (4 0) cosine coefficients.

[0079] If the magnet in use has a frequency of 7.692462 MHz and a spin echo sequence with a sampling frequency of 217 Khz, approximation processed by a polynomial like the one described above gives the following coefficients:

y=8.3385_(r) ⁴−25.7204_(r) ²

[0080] By applying the above, the following is obtained: $a_{20}^{\cos} = {{\frac{f_{sampling}\lbrack{MHz}\rbrack}{{f_{magnet}\lbrack{MHz}\rbrack}{2.56 \cdot 10^{- 4}}}\beta_{20}} = {\left. {\frac{0.00217}{7.682462{{\cdot 2}\text{,}5610^{- 4}}}\left( {{- 25}\text{,}7204} \right)}\Rightarrow a_{20}^{\cos} \right. = {{- 28}\quad {ppm}}}}$

[0081] In FIGS. 2A and 2B, in the XY plane, i.e. in the axial section plane, two baffles are positioned, to form an angle of 54.7360 with the Y axis. This suppresses the effect of distortion of the (2 0) cosine coefficient.

[0082]FIG. 6 shows a chart of the (2 2) to (2 0) coefficients ratio as a function of the angular position of the baffles 5 and 5′. The peak apparently corresponds to the angle of 54.736°.

[0083] Here again, in order to detect distortions and to perform symmetrization, for each point the symmetric points with respect to the origin are considered, which are four when two baffles are provided.

[0084] In order to detect distortion for (2 2) cosine, the distortions of each point and of its symmetric counterparts are summed up and divided by 4. Similarly to the above, the chart of a distortion curve, being a function of the distance r from the center of the image is obtained. This curve may also be further symmetrized with respect to the origin to suppress odd terms. From the basic relation indicated in the previous paragraph and identified by numeral (1), which describes the ratio between the symmetrized distortion of the image and the inhomogeneity of the magnetic field, the following is obtained: ${\Delta \quad {X\lbrack{pixels}\rbrack}} \cong {{1.81 \cdot 10^{- 4}}\frac{f_{sampling}\lbrack{MHz}\rbrack}{{f_{magnet}\lbrack{MHz}\rbrack}{2.56 \cdot 10^{- 4}}}\left( {{a_{22}^{\sin}r^{2}} + {O\left( r^{4} \right)}} \right)}$

[0085] Approximation with the following polynomial

y=β ₁ _(^(r)) ²+β₂ _(^(r)) ⁴

[0086] gives the (2 2) cosine coefficient, by equating quadratic terms.

[0087] In order to compute the (2 1) sine coefficient, the following operation has to be performed, where A, B, C, D indicate the distortions of points A, B, C, D, symmetrical to A. $\frac{\left( {A - D} \right) - \left( {B - C} \right)}{4}$

[0088] Hence, the chart of the distortion curve, being a function of the distance r from the center of the image is obtained: ${\Delta \quad {X\lbrack{pixels}\rbrack}} \cong {{1.92 \cdot 10^{- 4}}\frac{f_{sampling}\lbrack{MHz}\rbrack}{{f_{magnet}\lbrack{MHz}\rbrack}{2.56 \cdot 10^{- 4}}}\left( {{a_{21}^{\sin}r^{2}} + {O\left( r^{4} \right)}} \right)}$

[0089] Here again approximation is performed by the following polynomial:

y=β ₁ _(^(r)) ²+β₂ _(^(r)) ⁴

[0090] And by equating quadratic terms the (2 1) sine coefficient is obtained.

[0091] With reference to FIGS. 3A and 3B, the latter show the test element 6, 6′ for measurement of the (2 2) sine coefficient.

[0092] In this case, two baffles are positioned, as in the example of FIGS. 2A and 2B, in the XY plane, i.e. in the coronal section plane, but they are disposed along the bisectors of the Cartesian coordinates XY.

[0093] Symmetrization of distortions is performed as described for (2 1) sine coefficient determination. In this case, the law of distortion provides the equation: ${\Delta \quad {X\lbrack{pixels}\rbrack}} \cong {{1.92 \cdot 10^{- 4}}\frac{f_{sampling}\lbrack{MHz}\rbrack}{{f_{magnet}\lbrack{MHz}\rbrack}{2.56 \cdot 10^{- 4}}}\left( {{a_{21}^{\sin}r^{2}} + {O\left( r^{4} \right)}} \right)}$

[0094] As described above, approximation is performed by the following polynomial:

y=β ₁ _(^(r)) ²+β₂ _(^(r)) ⁴

[0095] And by equating quadratic terms the (2 2) sine coefficient is obtained.

[0096] The (2 1) cosine coefficient is measured substantially as described above. In this case, baffles 7, 7′, like those 6, 6′ of FIGS. 3A and 3B, even as regards angular arrangement are positioned in the sagittal section plane, i.e. in the XZ plane, as shown in FIGS. 4A and 4B. The operations are the same, and the following law for describing the distortion curve is obtained: ${\Delta \quad {X\lbrack{pixels}\rbrack}} \cong {{1.92 \cdot 10^{- 4}}\frac{f_{sampling}\lbrack{MHz}\rbrack}{{f_{magnet}\lbrack{MHz}\rbrack}{2.56 \cdot 10^{- 4}}}\left( {{a_{21}^{\sin}r^{2}} + {O\left( r^{4} \right)}} \right)}$

[0097] By obtaining, as described herein, the coefficients of the most relevant degrees of the polynomial expansion in spherical harmonics of the field, it is possible to redetermine, on the basis of nominal coefficients, the correcting magnetic charges as well as the positions thereof in the magnetic structure, whereby the magnetic field homogeneity may be restored within predetermined tolerances, to suppress the coarsest distortions.

[0098] In FIG. 5, the above described steps are indicated in box E, IS, whereas the following steps for determining the magnetic charges or the ferromagnetic correcting elements M1, M2, . . . , Mn and their position as well as their actual magnetic structure positioning 1, 2, 3 on the magnetic structure are schematically indicated in boxes D, PO.

[0099] The operations for computing and determining distortions may be performed outside the machine, or the software for automatically determining correction magnetic charges and their position from the comparison of data detected in the measuring session with stored nominal data, in relation to the different test elements 4, 5, 6, 7, as shown in the schematic flow chart of FIG. 5, may be loaded in the Nuclear Magnetic Resonance machine.

[0100]FIGS. 7 and 8 show a preferred embodiment of a test element for performing the measurement and correction procedures indicated and detailed herebefore.

[0101] This test element has all the individual elements 4, 5, 5′, 6, 6′, 7, 7′ permanently therein. These may be mounted in a fixed manner in the test element, but may be also possibly removable or a set of individual test elements may be provided, each with the baffles in one of the configurations described in FIGS. 1, 2, 3, 4.

[0102] The test element 10 consists of a cylindrical Plexiglas container 1, wherein the baffles 4, 5, 5′, 6, 6′ and 7, 7′. The baffles 6, 7 and 6′, 7′ are identical, with the test element, resting on the circular base, only having to be rotated through 90°. Obviously, as illustrated and denoted by 7, 6 and 7′, 6′, there may be also provided a set of baffles in a plane perpendicular to the one containing the first baffles 6, 7 and 6′, 7′, so that the test element does not have to be moved or oriented several times, and that a single first positioning may be needed to perform any further measurement. The baffles 4, 5, 6, 7, 5′, 6′, 7′ extend along axes which form angles indicated with reference to FIGS. 1, 2, 3, 4, all intersecting in the origin 0, which is provided in coincidence with the median point of the cylindrical test element 10.

[0103] The baffles 4, 5, 6, 7, 5′, 6′, 7′ consist of thin diaphragms or ribs separating apertures 110 formed in a plexiglas disk mounted inside the test element 10.

[0104] As for baffles 6, 7, 6′, 7′, these are formed in diametral walls of the test element, in a single diametral wall if the test element is to be rotated for performing different measurements, or in two perpendicular diametral walls. The two walls are denoted by 210 and 210′ in FIG. 8. The baffles consist of thin ribs separating apertures 310, 410.

[0105] In the upper portion, i.e base of the test element 10, an opening is provided which may be sealed by a plug 11, whereby the cylindrical sealed container-like test element is filled with a liquid which generates Nuclear Magnetic Resonance signals, e.g. simple water , and which completely fills the test element, including the different apertures 110, 310, 410 which define the baffles 4, 5, 5′, 6, 6′, 7, 7′.

[0106] Upon measuring, the test member is inserted inside the receiving coil of a Nuclear Magnetic Resonance imaging machine in such a position that the axis of the cylindrical member is oriented perpendicular to the static magnetic field indicated by the arrow Bz in FIG. 8. The member is positioned with the baffles, i.e. walls 210, oriented parallel to the two Cartesian axes which define the planes perpendicular to the direction of the static magnetic field, i.e. the faces of the two opposite poles 2, 3, whereas the baffle 4 is oriented along the y axis of said plane, i.e. according to one of the reading gradients.

[0107] The member has extensions 510, 610, 710 for easy positioning and grasping thereof.

[0108] The member as shown in FIGS. 7 and 8 is particularly suitable for measurement of field coefficients as described hereinbefore with reference to the method and to FIGS. 1 to 5.

[0109] Obviously, in combination with the test member, a Nuclear Magnetic Resonance imaging machine has to be provided which has at least one magnetic structure for generating static magnetic fields with a predetermined orientation; electronic means and coils for generating spin echo sequences, as well as selecting, reading and encoding gradients; coils for receiving echoes and an electronic signal processing unit, said electronic unit being programmed for performing at least some of the steps of the method described hereinbefore, and there being provided, in combination therewith, at least one set of correcting magnetic charges, whose values are fractions or multiples of magnetic charge units.

[0110] The invention is not limited to what has been described and illustrated herein, but may be widely varied, especially as regards construction, without departure from the guiding principle disclosed above and claimed below. 

1. A method for performing verification measurements and corrections of magnetic fields in magnets for Nuclear Magnetic Resonance imaging machines, providing the following steps: making a test member (10; 4, 5, 5′, 6, 6′, 7, 7′) having such characteristics as to generate predetermined and known Nuclear Magnetic Resonance images (I); obtaining a Magnetic Resonance image from said phantom (I′); comparing the known theoretical image with the detected image (I+I′) and/or the received signals wherefrom the detected image is reconstructed with the corresponding set of signals, related to the known theoretical image in any of their processing or imaging steps; determining the differences between the detected image and the known image (I+I′, IS) or between the set of detected signals and the set of signals corresponding to the known image; computing correction parameters, i.e. the correcting magnetic charge (M1, M2, . . . , Mn) and/or the number of correcting magnetic charges having a predetermined value, or the volume of ferromagnetic material and/or the number of ferromagnetic correcting elements of a predetermined volume, and the position of said correcting elements on the magnetic structure (1). Characterized by the further steps: Choosing a mathematical theoretical function describing the field, particularly a polynomial expansion; Providing a phantom with a structure suited to the said function describing the filed and which produces one or more selected images, each of these images being correlated to one or to a limited number of selected low order coefficients of the function describing the filed, such that the differencies between a theoretical image of the phantom and the real image of the phantom depends only to the one or to the limited number of selected coefficients of the function describing the magnetic field.
 2. A method as claimed in claim 1 , characterized in that the defined mathematical structure for describing the field is choosen as a representation of the field in spherical harmonics, and such that only different orders and degrees of harmonics and hence of coefficients of the field may be examined separately.
 3. A method as claimed in claim 1 or 2 , characterized in that it provides that the test member (10) has elements (4, 5, 5′, 6, 6′, 7, 7′) which do not emit Nuclear Magnetic Resonance signals and are related to a definite harmonic or to precise coefficients of a certain order of a mathematical description of the field corresponding to a polynomial expansion, particularly by spherical harmonics, there being provided means for application of so-called reading gradients of the magnetic field, which only detect the echo signal along certain directions, said directions being selected in such a manner as to suppress the contributions from the magnetic field described by coefficients other than the ones being examined, to isolate the contribution of the coefficients being examined from the other remaining coefficients.
 4. A method as claimed in one or more of the preceding claims, characterized in that the describing polynomial expansion is selected in spherical harmonics.
 5. A method as claimed in one or more of the preceding claims, characterized by the following steps: a) detecting the image (I′) of the predetermined test element (10, 4); b) symmetrizing deviations between the actually detected image (I′) and the nominal ideal one (I) with respect to the center of the image and/or to the predetermined origin of the coordinate system for the mathematical description of the magnetic field; c) defining a curve of symmetrized deviations having the relevant field coefficients as a variable; d) determining a polynomial for approximating the curve of symmetrized deviations, deriving from the differences between the actual image (I′) and the ideal, theoretical image (I); e) determining the coefficients based on the system of the two consistent equations of steps c) and d); f) computing, from the mathematical description in spherical harmonics of the field, such number, magnitude and position on the magnetic structure (1) of the correcting elements as magnetic charges or volumes of ferromagnetic material (M1, M2, . . . , Mn), as to bring the measured values of the field coefficients back within the nominal values corresponding to the homogeneity characteristics of the field required for detecting useful Nuclear Magnetic Resonance images; g) manually positioning such correcting elements (M1, M2, . . . , Mn) on the magnetic structure (1, 2, 3).
 6. A method as claimed in one or more of the preceding claims, characterized in that, once correction has been performed, steps a) to f) may be repeated to verify the effectiveness of said correction and, when necessary, an additional correction step g) may be performed.
 7. A method as claimed in one or more of the preceding claims, characterized in that the elements (4, 5, 5′, 6, 6′, 7, 7′) for verifying the different coefficients may be progressively mounted on the test member (10) or fixedly provided thereon.
 8. A method as claimed in one or more of the preceding claims, characterized in that each of these test elements (4, 5, 5′, 6, 6′, 7, 7′) corresponds to at least one harmonic or to at least one set of coefficients and are made consistently with the selected mathematical description of the field and/or eventually to the symmetries of the magnetic structure (1).
 9. A method as claimed in one or more of the preceding claims, characterized in that it provides a manual graphic comparison between the detected images (I′) and the theoretical images (I) wherefrom the mathematical computation is performed to determine correction data (D).
 10. A method as claimed in one or more of the preceding claims, characterized in that the machine has a memory wherein the theoretical comparison image (I) is loaded.
 11. A method as claimed in one or more of the preceding claims 1 to 9 , characterized in that it provides that the machine determines quantitatively the differences between the ideal image (I) and the detected image (I′) by computation of the distance in pixels of the point or portion or area of the detected image (I′) from the position of the corresponding point or area of the ideal image (I).
 12. A method as claimed in one or more of the preceding claims, characterized in that it provides the use of a test element (10) whose undetectable portion has the form of a rib, step or baffle (4) having a longitudinal orientation parallel to the static magnetic field (Bz), i.e. disposed in the axial section plane of the magnetic field, when the latter is described by spherical harmonics, said baffle being related to the detection of (2 0) cosine and (4 0) cosine coefficients only, with reference to a mathematical description in spherical harmonics.
 13. A method as claimed in claim 12 , characterized in that it provides the application of so-called reading field gradients oriented perpendicular to the baffle.
 14. A method as claimed in one or more of the preceding claims, characterized in that two crossed baffles (5, 5′) are used as a test element, which form an angle of 109.472°, symmetrically with respect to one of the two axes describing the axial section plane of the magnetic field (Bz), which baffles are related to measurement of contribution of the field (Bz) only relating to (2 2) cosine and (2 1) sine coefficients.
 15. A method as claimed in one or more of the preceding claims, characterized in that two crossed baffles (6, 6′) are used as a test element, forming an angle of 90°, symmetrically with respect to an axis which describes the coronal section plane of the static magnetic field (Bz), and being related to measurement of (2 2) sine coefficients only, whereas, by the same configuration of baffles (7, 7′), though disposed in the sagittal section plane, field contributions of (2 1) cosine coefficients are only detected.
 16. A test element for verification measurements and corrections of magnetic fields, particularly in Nuclear Magnetic Resonance machines, characterized in that it comprises, separately or in combination, at least one of the baffles (4, 5, 5′, 6, 6′, 7, 7′) as claimed in one or more of the preceding claims 12, 13, 14, 15, or all aforesaid baffles at the same time.
 17. A test element as claimed in claim 16 , characterized in that said baffles (4, 5, 5′, 6, 6′, 7, 7′) are made of a plastic material, particularly of the so called Plexiglas.
 18. An element as claimed in claim 16 or 17 , characterized in that it includes a central rectilinear baffle (4) which is crossed, with reference to the plane parallel to a peripheral edge, by two crossed baffles (5, 5′), passing through the median area of the rectilinear baffle (4) and forming, symmetrically with respect to the rectilinear baffle (4), an angle of 109.472°, whereas two additional crossed baffles (6, 6′) are provided in the plane perpendicular to the rectilinear baffle (4) and containing the other sides thereof, which form, symmetrically with respect to said rectilinear baffle (4), an angle of 90° and intersect the rectilinear baffle in the central area, yet two additional baffles (7, 7′), crossed at 90° being provided in another plane perpendicular to the plane containing the two baffles crossed at 90°.
 19. An element as claimed in one or more of the preceding claims 16 to 18 , characterized in that the baffles are formed by thin walls separating apertures (110, 310, 410) in axial or diametral walls (210) intersecting along the three main directions of a Cartesian system.
 20. An element as claimed in one or more of the preceding claims, characterized in that the various baffles (4, 5, 5′, 6, 6′, 7, 7′) and/or the walls (210) wherein they are formed are mounted permanently in a plastic, particularly Plexiglas, container (10), which is filled with a liquid or a substance emitting a Nuclear Magnetic Resonance signal, for instance water, the container being water-tight and provided with a sealable filling opening (11).
 21. An element as claimed in one or more of the preceding claims 16 to 20 , characterized in that it comprises a water-tight cylindrical container (10).
 22. An element as claimed in one or more of the preceding claims, characterized in that the baffle/s are made of a material suitable for generating Nuclear Magnetic Resonance signals, whereas the filling material is such that it does not emit Nuclear Magnetic Resonance signals, the generation of separator surfaces of a predetermined shape, for separating contrasting image areas, being relevant.
 23. A Nuclear Magnetic Resonance imaging machine, comprising a magnetic structure (1) for generating static magnetic fields, which defines a cavity for accommodating bodies or body parts to be examined; a coil for generating sequences for exciting, selecting and reading Nuclear Magnetic Resonance echoes; a coil for receiving said echoes; a unit for processing the echo signals received, to generate digital images therefrom to be displayed on display means, characterized in that the processing unit is of the programmable type and has at least one memory which stores an algorithm for measuring field coefficients according to the method as claimed in one or more of the preceding claims 1 to 15 , whereas, in combination with said machine, there is provided a test element as claimed in one or more of the preceding claims 16 to 21 , whose ideal image may be stored in the processing unit of the machine, and at least one set of correcting magnetic charges (M1, M2, . . . , Mn). 